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ABSTRACT 

We present 2D MHD simulations of the radiatively driven outflow from a rotating hot star with 
a dipole magnetic field aligned with the star's rotation axis. We focus primarily on a model with 
moderately rapid rotation (half the critical value), and also a large magnetic confinement parameter, 
rj* = B^Rl/MVoc = 600. The magnetic field channels and torques the wind outflow into an equatorial, 
rigidly rotating disk extending from near the Kepler corotation radius outwards. Even with fine-tuning 
at lower magnetic confinement, none of the MHD models produce a stable Keplerian disk. Instead, 
material below the Kepler radius falls back on to the stellar surface, while the strong centrifugal force 
on material beyond the corotation escape radius stretches the magnetic loops outwards, leading to 
episodic breakout of mass when the field reconnects. The associated dissipation of magnetic energy 
heats material to temperatures of nearly 10 s K, high enough to emit hard (several keV) X-rays. Such 
centrifugal mass ejection represents a novel mechanism for driving magnetic reconnection, and seems 
a very promising basis for modeling X-ray flares recently observed in rotating magnetic Bp stars like 
a Ori E. 

Subject headings: MHD — stars: winds — stars: magnetic fields — stars: rotation — stars: flare — 
X-rays: stars 



1. INTRODUCTION 

Magnetic fields can greatly influence stellar winds by 
guiding the outflowing plasma along fi eld lines (e.g., 
iMacGregorl Il988t iShore fc Brownl 11990(1 . The extent 
of this influence depends largely on the magnetic field 
strength and the wind mass loss rate. Two-dimensional 
(2D) magnetohydrodynamical (MHD) simulations of ax- 
isymmetric hot-star winds with no rotation show that the 
competition between the wind and the field can be quite 
well characterized by a singl e wind magnetic confinemen t 
parameter, 77*, as defined bv lud-Doula fc Owocki {2002), 
and described further in jJS] below. 

For significant confinement (i.e., 77* > 10), mate- 
rial driven from opposite footpoints of closed magnetic 
loops is forced to collide near the loop tops, leading 
to a Magnetically Confined Wind Shock (MCWS) sce- 
nario for production of inte rmediate-hardness X-rays 
l|Babel fc MontmerTel U997aTbT) . For the magnetic hot 
star 9 1 Ori C, which is inferred to have both a mod- 
erately strong wind (with mass loss rate « 10 -7 M Q /yr) 
and a large-scale, d ipole field of strength as 1100 G 
l|Donati et alj |2002l), recent MHD simulations based on 
this MCWS paradigm provide quite good ag reement with 
the X -ray spectrum observed by Chandra l|Gagne et al.l 
2005). Indeed, because 9 1 Ori C is a relatively slow rota- 
tor (with period ~ 15 days), even the observed rotational 
modulation can be well accounted for by simply varying 
the observational perspective within a non-rotating dy- 
namical simulation model. 

The purpose of this letter is to report on initial 
MHD simulation results for more rapidly rotating hot 
stars, in which the rotation can have a direct dynam- 
ical effect on the overall evolution of the wind outflow 
and its confinement. The results have particular rele- 
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vance in the context of two proposed models for rotat- 
ing, magnetic hot stars, known genera lly as the Mag- 
netic ally Torqued Disk (MTD) model l|Cassinelli et alJ 
2002) and the Rigidly Rotating M agnetosphere (RRM) 
model ijTownsend fc Owockill200l . The MTD has been 
proposed as a scenario for producing Keplerian decre- 
tion disks, such as inferred from the Balmer line emis- 
sion of classical Be stars. The RRM leads instead to 
rigid-body disks or clouds that provide a promising way 
to explain the rotationally modulated Balmer emission 
from str ongly magnetic Ap/Bp stars like the B2pV star 
a Ori E <|Townsend et al.ll2005[l . 

After first casting the relative importance of magnetic 
confinement vs. rotation in terms of characteristic Alfven 
vs. Kepler corotation radii (©, we show in ^Ithat even 
for parameters fine-tuned to optimize the chances of a 
producing a Keplerian disk, our 2D MHD simulations 
for rotationally aligned, dipole fields instead produce a 
combination of inflow and outflow for material below and 
above the Kepler corotation radius, and never lead to the 
kind of extended Keplerian disk envisioned in the MTD 
scenario. 

On the other hand, we show in Sj3|that a model with 
strong magnetic confinement (77* = 600) does lead to a 
rigid-body disk, much as predicted in the RRM model. 
However, the outer regions of this disk are character- 
ized by episodic breakout of material, driven by the 
strong net outward centrifugal force associated with the 
spin-up of material into corotation with the underlying 
star. Moreover, the energy release associated with mag- 
netic reconnection in these breakout events heats mate- 
rial to extremely high temperatures, providing a promis- 
ing mechanism for explaining the very hard X-ray flare 
emission recently observed from magnetic Bp stars such 
as a Ori E. We conclude in jj3 with a brief summary 
of our findings, noting also the need for further work on 
the detailed physics of this somewhat-novel, centrifugally 
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driven reconnection mechanism. 

2. ALFVEN VS. KEPLER COROTATION RADII 
MHD simulations (jud-Doula fc Owockil 120021: 



lud-Doulal I2003D indicate that the effectiveness of 
magnetic fields in channeling a stellar wind outflow can 
be characterized by the ratio between the magnetic and 
wind energy densities 
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This ratio furnishes an indication of whether the mag- 
netic field at radius r dominates the wind outflow (77 3> 
1), or vice versa (77 <C 1). As indicated by the second 
equality, its overall scaling is set by the "magnetic con- 
finement parameter" 77* = B 2 R 2 /MVqo, which depends 
on the strength of the field B* at the stellar surface ra- 
dius r = i?* , and on the wind terminal momentum MVoo . 
The radial variation of 77 is then modeled in terms of a 
magnetic power-law index q (= 3 for a dipole) and a ve- 
locity index /3 (« 1 for a standard CAK wind). If for 
simplicity we ignore the wind velocity variation (i.e. by 
taking f3 — 0), the Alfven radius - at which the magnetic 
and wind energy densities are equal, t](Ra) = 1 — can be 
solved explicitly as 



Ra = vl /4 R* 



(2) 



As shown by simulati on results summarized in 
lud-Doula fc Owockil l)2002T) , this Alfven radius provides a 
reasonable estimate for the maximum radius of magnetic 
loops that remain closed in a wind outflow. Moreover, 
since in rotating models such closed loops tend to keep 
the outflow in rigid-body rotation with the underlying 
star, it also defines the radius of maximum rotational 
spin- up of the wind azimuthal velocity v$. 

To characterize such rotational effects, let us next de- 
fine a Kepler corotation radius Rk- Analogous to the 
terrestrial geostationary orbital radius, Rk is the point 
at which the centrifugal force arising due to corotation 
balances the local gravitational force from the star. That 
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where V eq is the equatorial rotation velocity of the star. 
Solving this equation, we find the Kepler radius as 



R K = W' 
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where W = V cq /V cr it is ratio of the star's equatorial ve- 
locity to the critical velocity V cr it = yjGM*/R*. 

Finally, it is also worth noting here that corotation out 
to an only slightly higher escape radius, 



R E = 2 1/3 R K = 2 1/3 W 
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would imply an azimuthal velocity that equals the 
local escape speed from the star's gravitational field. 

3. ATTEMPTS TO FINE TUNE FOR A KEPLERIAN DISK 

The above scalings suggest that a likely necessary con- 
dition for propelling outflowing material into a Keple- 
rian disk (as envisioned in the MTD scenario for classi- 
cal Be stars) is to choose a combination of parameters 
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Fig. 1. — The key domains in a parameter plane of magnetic 
field strength (as represented by yJr\Z ~ B„) vs. stellar rotation 
(as represented by the critical rotation fraction W ~ V cq ). The 
black area to the lower left denotes models in which the magnetic 
spin-up of the azimuthal velocity is insufficient to reach the 
velocity required for Keplerian orbit, while the white area in 
the upper right denotes models in which the spin-up leads to 
that exceed the escape speed i>e- Only in the gray area, bounded 
between curves for models defined by Rj± = _Rk an d Ra = Re, 
does spin-up into Keplerian orbit seem possible. The dot positioned 
in the middle of this area indicates the parameters of the MHD 
simulations described in J3] 



for magnetic confinement vs. stellar rotation such that 
-Rk < Ra < Re- In the parameter plane of ^/fj^ ~ £?* 
vs. W ~ V eq illustrated in Fig. ^ the required combina- 
tion is represented by the relatively narrow gray band. 

To obtain a model that has the best chance of yielding 
a Keplerian disk, we therefore conduct 2D axisymmetric 
MHD simulations of a star having parameters 77* = 10 
and W — 0.5 that fall in the middle of the gray area in 
Fig- ID The initial magnetic field geometry is assumed to 
be a dipole aligned with the rotation axis. Calculations 
are u ndertaken using the p ublicly available ZEUS-3D 
code l)Stone fc Normanll 1992(1 . augm ented to include ra - 
diative driving and radiative cooling ijOagne et all 2005) , 
with typical early B-star values for the stellar radius and 
mass, i?* = 6.4 i?0 and = 11 M®. The numeri- 
cal grid and boundary con ditions follow those given in 
lud-Doula fc Owockil |2b02). Since the star is rotating 
at one half its critical value, the oblateness and gravity 
darkening of the stellar surface are relatively small, and 
are neglected in our simulations. 

Figure [21 illustrates results from the simulation. The 
left-hand panel shows the logarithm of density (super- 
imposed with field lines) at a time 90 ks after the intro- 
duction of the dipole field into a spherically symmetric 
wind. Superficially, these conditions do temporarily re- 
semble the MTD scenario. However, a closer inspection 
shows that most of the density-enhanced equatorial re- 
gion does not have the appropriate necessary for a 
stable, stationary, Keplerian orbit. Indeed, in just a few 
ks of subsequent evolution this putative "disk" becomes 
completely disrupted, characterized generally by infall of 
the material in the inner region, i.e. below the Kepler 
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Fig. 2. — Density in a 2D MHD simulation with rj, = 10 and 
W = 1/2, shown at time snapshots of 90 ks (left) and 390 ks (right) 
after a dipole field is introduced into an initially steady, spherical, 
unmagnetized, line-driven stellar wind. The curves denote mag- 
netic field lines, and the vertical dashed lines indicate the equatorial 
location of the Kepler (%), Alfven (-Ra)i and escape (-Re) radii 
defined in eqns. J1J, J5Ji an d O > respectively The arrows denote 
the upward and downward flow above and below the Kepler radius, 
emphasizing that the material never forms a stable, orbiting disk. 



radius Rk, and by outflow in the outer region above the 
Kepler radius. The right-hand panel demonstrates this 
by showing the irregular form of the dense compression 
at an arbitrarily chosen later time (390 ks). The bold 
arrows emphasize the flow divergence of the dense mate- 
rial both downward and upward from the Kepler radius. 
(The overall evolution is most vividly illustrated through 
animations available on the internet 2 .) 

We have carried out a number of other similar MHD 
simulations exploring parameter space in both the rota- 
tion rate and magnetic field strength directions. In every 
case, we find that the equatorial compressions are dom- 
inated by similar radial inflows and/or outflows, with 
no stable Keplerian disks forming. We thus conclude 
that stationary, large-scale magnetic fields are dynami- 
cally ill-suited for producing Keplerian or near-Keplerian 
decretion disks. 

4. A CENTRIFUGAL BREAKOUT MECHANISM FOR 
X-RAY FLARING 



A recent study by iTownsend fc Owockil (j2005f) sug- 
gests that, in the limit of very strong magnetic con- 
finement (77* 3> 1), magnetic torquing can readily pro- 
duce rigid-body disks. In this strong- field limit, the field 
lines can be idealized as rigid pipes that channel the 
wind outflow and force it to corotate with the underly- 
ing star. This leads to a rigidly rotating magnetosphcrc 
in which material accumulates near minima of the ef- 
fective gravito-c e ntrifu gal potential along each field line. 
ITownsend et alJ l|2005|) show that the resulting circum- 
stellar clouds provide a natural explanation for the ob- 
served rotational modulation of both the photometric 
continuum and Balmer line emission in the B2pV star 
a Ori E. 

A key feature of this RRM model is that the magnetic 
field not only torques up outflowing material into corota- 
tion around the star, but then also holds it down against 
the net outward centrifugal acceleration beyond the Ke- 
pler radius. In the semi-analytic RRM analysis, the field 
is idealized as being strong enough to enforce corotation 
and confinement at arbitrarily large radii, and for arbi- 
trarily long times. More realistically, however, for large 
but finite field strengths, the secular accumulation of ma- 
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terial into the RRM means that these idealizations must 
break down once the total outward centrifugal force ex- 
ceeds a level that can be balanced by the finite magnetic 
t ension. 

ITownsend fc Owockil l)2005l) speculate that this will 
lead to centrifugal breakout of material, and give esti- 
mates for the expected characteristic breakout timescale 
as a function of distance above the Kepler radius (see 
their Appendices A and B). For the case of a Ori E, 
they estimate that the largest-scale evacuations can be 
expected every ~ 10 — 100 yr, but they also anticipate 
a whole hierarchy of breakout events extending down to 
timescales of days. 

Our present MHD simulations allow us to examine the 
dynamics of this breakout process. Unfortunately, at the 
confinement parameter 77* ~ 10 7 appropriate to a Ori E, 
the Alfven speed becomes very large, implying a very 
small Courant time step to ensure numerical stability; 
this together with the full 3D geometry needed to model 
the oblique-field tilt relative to the rotation axis makes 
direct MHD simulation of the a Ori E case prohibitively 
expensive. Instead, we consider the more tractable, 2D 
axisymmetric case of a field-aligned dipole with only 
moderately strong magnetic confinement (77* = 600). As 
appropriate for a Ori E, we again set the rotation to half 
the critical rate (W = 0.5), with other stellar and grid 
parameters also the same as used in S|21 

Snapshots from the simulations, plotted in Fig. El re- 
veal that there forms a rigidly rotating disk near and 
above the Kepler radius Rk/R* ~ 0.5~ 2 / 3 rs 1.6, as pre- 
dicted by the RRM model. But in the outer regions of 
this disk (near the Alfven radius Ra/R* ~ 600 1//4 ~ 5), 
there is indeed also a semi-regular sequence of break- 
out events. During a breakout, the magnetic field lines 
become so drawn out by the ejected plasma that they 
reconnect, snapping back toward the star (cf. middle 
and right panels of Fig. [3J| . The energy release associ- 
ated with this reconnection, and its subsequent dissipa- 
tion via radiative cooling, represents a strong candidate 
for producin g the hard X-ray fl ares observed in a Ori E 
bv ROSAT fcroote fc Schmittll200^ and XMM-Newton 
l|Sanz-Forcada et al.ll200"l ~ 

This scenario is reminiscent of fast magne tic recon- 
nection processes invoked for solar flares (e.g.. iPetschekl 
1966§ Eitti|2003). But in the present case, the work to 
stretch the magnetic fields comes from the rotational en- 
ergy of the dense plasma at the edge of the disk, with 
the trigger for reconnection being the wind impingent 
at an oblique angle from higher latitudes. The centrifu- 
gal energy stored temporarily in the magnetic field is 
thereby converted into thermal energy, heating nearby 
material to temperatures T ~ 10 s K, high enough to pro- 
duce the hard (> 2 keV) components of the X-ray flares 
observed in a Ori E. A similar mechanism may also be 
relevant to X-ray outbursts recently observed in other 
(possibly magneti c) early-type stars, su ch as the mid-B 
star HD 38563S ( Yanagida et all 120041) and the late-0 
star 2 Ori A l)Feigelson et al.ll2002|) 

The centrifugal breakout here provides an inter- 
esting c ounterexample to the "slingshot prominence" 
model of Uardine & van Ballegooiienl l)2005[) . This cool- 
star, coronal mode l envisions a tearing- mode instability 
llFurth et alJ 11963) in the current sheet dividing open 
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Fig. 3. — MHD simulation results for a radiativcly driven stellar wind from a star rotating at half the critical rate, with a rotation- 
aligned dipole field that has a magnetic confinement parameter r?» = 600. The grayscales show the logarithm of density (top panels) 
and temperature (bottom panels) at three sample time snapshots, chosen to show the centrifugal ejection of material accumulated in an 
equatorial disk. The associated stretching and eventual reconnection of the field lines heats material to temperatures of order 10 8 K, hot 
enough to emit hard (few keV) X-rays. 



field regions in an overlying, gas-pressure-driven wind. 
This leads to a magnetic reconnection that effectively 
traps some of the outflowing wind material, contribut- 
ing then to an exten ded, quasi-static coronal prominence 
ijJardine et alJl200l|) . Although the centrifugal support 
of this prominence is similar to our RRM scenario for hot- 
stars, in this cool-star case the high coronal back-pressure 
chokes off any footpoint upflow that could continuously 
feed the accumulation region above the Kepler radius. 
In our hot-star model, it is this secular accumulation of 
material that eventually forces centrifugal breakout and 
drives the associated reconnection. As such, unlike the 
tearing-mode reconnection postulated in this cool-star 
context, we believe the overall rate of reconnection and 
mass ejection should be relatively independent of the spe- 
cific instabilities and plasma processes at play. 

In real plasmas, the details for magnetic reconnec- 
tion can be quit e complex, involving non-thermal parti- 
cle acceleration dWiechen et al.lll995fl and /or high -speed 
jets that form dissipative shocks l)Petscheklll966|) . with 
the distribution of heat modifie d by thermal conduction 
or no n-local particle transport (Mivagos hi fc Yokovamal 
2003). The Zeus-3D numerical code used here is based 
on ideal an MHD formalism that effectively ignores these 
complexities, and indeed does not even include any ex- 
plicit resistivity that could heat material via Ohmic dis- 
sipation. Instead, magnetic reconnection is enabled by 
the finite spatial grid resolution, leading to grid-scale 
currents inducing high-speed jets that quickly dissipate 
the flow energy in localized shock fronts. Nonetheless, 
the reconnection rates observed in o ur models are con - 
sistent with kinetic simulations (e.g.. lShav et al.l fl999). 
and the code still approximately fulfills global energy 



conservation, thus providing a first description of the 
net dissipation of magnetic energy into plasma heating. 
Hence, we believe the simulations presented here repre- 
sent a physically reasonable scenario for heating by cen- 
trifugally driven reconnection, with the associated X-ray 
emission providing a promising explanation for the X-ray 
flares observed from magnetic Ap/Bp stars. 

5. SUMMARY 

The rigid-body corotation makes the equatorial disks 
in our simulations quite distinct from the Keplerian 
disk envisioned in t he MTD scenario promoted by 
ICassinelli et al.l (J2002). We find no evidence of an ex- 
tended region of Keplerian rotation. Instead the re- 
sults correspond more closely to thos e deri ved in the 
RRM analysis of iTownsend fc Owockil (|2005Tl . although 
the confinement parameter we consider is near the lower 
end of the range of validity for the basic rigid-field ap- 
proach. 

Moreover, our numerical simulations further suggest a 
novel mechanism for magnetic reconnection, driven by 
episodic centrifugal breakout events. The strong recon- 
nection heating associated with such breakout represents 
a promising mechanism for explaining X-ray flares ob- 
served in magnetic Bp stars such as a Ori E. 

In future work, we plan to use more-extensive numer- 
ical tools to investigate in greater detail the reconnec- 
tion processes within this centrifugally driven breakout 
scenario. We also plan to explore hybrid methods that 
could allow extensions of MHD simulations to larger con- 
finement parameters, including 3D models of tilted-field 
cases more appropriate to a Ori E and other magnetic 
hot-stars. 
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